Fluctuation— dissipation relations in the activated regime of simple strong— glass models 
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We study the out-of-equilibrium fluctuation-dissipation (FD) relations in the low temperature, 
finite time, physical aging regime of two simple models with strong glass behaviour, the Fredrickson- 
Andersen model and the square-plaquette interaction model. We explicitly show the existence of 
unique, waiting-time independent dynamical FD relations. While in the Fredrickson-Andersen 
model the FD theorem is obeyed at all times, the plaquette model displays piecewise linear FD 
relations, similar to what is found in disordered mean-field models and in simulations of supercooled 
liquids, and despite the fact that its static properties are trivial. We discuss the wider implications 
of these results. 
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The common feature to all glassy systems, like super- 
cooled liquids, spin glasses, and, to a certain degree, even 
soft materials like gently driven powders, is an extremely 
slow relaxational dynamics at low temperatures or high 
densities (for reviews see Often, glassy systems 

display a dynamical behaviour known as aging, which 
corresponds to the asymptotic regime in which one-time 
quantities, like energy or magnetization, are stationary, 
but two-time quantities, like response and autocorrela- 
tion functions, still depend on the time elapsed since 
the system was prepared, rather than just on time dif- 
ferences, as in equilibrium. While in this situation re- 
sponse and autocorrelations do not obey the fluctuation- 
dissipation theorem (FDT) exact results for mean- 
field models [^,|| have suggested that FD relations are 
generalized in a well defined way, that the breakdown of 
FDT can be understood in terms of 'effective' tempera- 
tures 1^, and that there may be a close connection be- 
tween out-of-equilibrium FD relations and equilibrium 
properties Nontrivial asymptotic FD relations have 
also been found in other (non-glassy) out-of-equilibrium 
situations, like ferromagnetic domain growth Q and sys- 
tems at criticality l^, and FD plots similar to the ones 
of discontinuous mean-field models have been observed 
in simulations of supercooled liquids pO| , pl] | and of frus- 
trated 10 and constrained lattice gases. 

However, both in experiments and simulations, a rel- 
evant regime is that of long but finite times where one- 
time quantities are not stationary, but are slowly relax- 
ing towards their equilibrium values, a situation known as 
'physical' aging [|l^ . A second issue is that near the glass 
transition activated processes, which are explicitly ex- 
cluded in mean-field, play an essential role. With respect 
to this, simulations of simple models with activated dy- 
namics JT5| , p^ have shown nonmonotonic response func- 
tions, which, superficially, lead to meaningless FD rela- 
tions (something analogous occurs in models of vibrated 
granular matter ||l^). And finally, there is evidence, at 
least for molecular glasses, for the absence of any ther- 
modynamic phase transition underlying the dynamical 
arrest [^, so that if nontrivial FD relations exist for 
these systems they cannot be interpreted in terms of the 



structure of an equilibrium glass phase. 

The purpose of this Letter is to address the problem 
of whether well defined out-of-equilibrium FD relations 
can be obtained when all of the above factors are taken 
into account. We do this by considering the case of two 
simple non-disordered or frustrated systems with triv- 
ial statical properties but dynamical behaviour charac- 
teristic of a strong-glass [0 (like Arrhenius relaxation 
time, exponential relaxation functions, etc.): the ID and 
2D Fredrickson-Andersen model [p^| j20|Jl5| ,|2l|| and the 2D 
square-plaquette interaction model [p2[. We show that 
for an appropriately defined class of observables non- 
trivial, unique out-of-equilibrium FD relations exist in 
the physical aging regime of these models. Using simple 
scaling arguments we find that in the FA model FDT is 
obeyed at all times. Remarkably, in the plaquette model 
the existence of two relevant timescales for the relaxation 
leads to piecewise linear FD relations, despite its trivial 
statics. 

Let us consider first the case of the Fredrickson- 
Andersen (FA) model |l9|, which corresponds to Ising 
spins = ±1 (i = l,...,iV), with Hamiltonian H = 
rii, where rii = (l-f(Ti)/2, and subject to asingle spin- 
fiip dynamics with the kinetic constraint that only spins 
which have at least one nearest neighbour in the up state 
are allowed to flip. This dynamics obeys detailed bal- 
ance, and the equilibrium is trivial. The energy density 
is given by the concentration c of up spins (or 'defects'), 
which in equilibrium becomes Coq = 1/(1 + e^^'^). At low 
temperatures c is very small, and since defects facilitate 
the dynamics, the system slows down. Isolated defects 
are locally stable and the system has to overcome energy 
barriers to evolve. There is a single activation barrier to 
the diffusion of defects AE = 1, which implies that re- 
laxation times follow the Arrhenius law tfa gConst./T^ 
characteristic of strong glass behaviour |2l[] . We wish to 
study the dynamical FD relations for long but finite times 
after a quench from Tq = oo to a low temperature T. 
This is the physical aging regime of the system in which 
one-time quantities change slowly with time as they relax 
towards their equilibrium values. For low temperatures, 
the concentration of defects c{t) = N^^ devel- 
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ops a plateau at Cp, which becomes longer the lower the 
temperature, and corresponds to the onset of activation. 
The plateau is reached at t 1 and this initial transient 
is T independent. We are interested in times i > 1, that 
is, from the plateau onwards, for which the relaxation 
proceeds through activated processes. 

The integrated response for an individual spin can be 
measured by the now standard method |^ of applying 
a random field perturbation, SH{t) = — e^ni, 
where Si — ±1 are iid random variables. For a con- 
stant field applied from waiting time w onwards, the 
integrated response at time t is given by x(t, w) = 
{Nh)~^ J2i^i{''^ii^))h, where {■)h indicates dynamical av- 
erage in the presence of the perturbation, and (•), av- 
erage over the realization of the random field. In Fig. 
1 (top left) we show w) for low temperatures and 
various waiting times for the ID FA model |^^, the be- 
haviour of the 2D version being similar. Notice that the 
response function is non-monotonic . The peak in the 
response decreases with increasing waiting time, and the 
function eventually becomes monotonic when the system 
equilibrates. On the other hand, the spin autocorrela- 
tion function C{t, w) = N^^ 'Y^^{'ni{t)ni{w)) is monoton- 
ically decreasing, as expected. The combination of non- 
monotonic response and monotonic correlation was in- 
terpreted as leading to meaningless FD relations for this 
model 15 1 . A more careful analysis reveals something 
different. Only defects contribute to the response func- 
tion which implies that w) should be proportional to 
c(t). Moreover, to leading order in e~^/'^ , a spin which 
can respond is isolated and, thus, in local equilibrium. 
This leads to the following scaling form for the response 
for T < 1 and m; > 1, 



X[t,W) Xoq 



(t-w) 



= 1/T 



(1) 



where Xoq is the equilibrium response function, which 
depends only on the time difference, and scales with the 
temperature through the relaxation time of up spins. Fig. 
1 (top right) presents the excellent collapse of all the data 
of Fig. 1 (top left) under Eq. (0). The above expression 
also reveals the nature of the non-monotonicity in the 
out-of-equilibrium susceptibility: it is the product of a 
decreasing function, c, and an increasing one, Xeq- Simi- 
lar arguments give the scaling form of the autocorrelation 
C{t,w): of the up spins at the earlier time, c{t) will re- 
main and c{w) — c(t) will disappear. Assuming that the 
autocorrelation of the former is proportional to the equi- 
librium one, and that they are uncorrelated to the latter, 
we obtain. 
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c{t)[c{w) - c{t)], (2) 



where Ccq is the equilibrium autocorrelation function. 
Fig. 1 (bottom left) shows the excellent collapse of auto- 
correlation functions under (^. 
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Top left: integrated response xi't, ™) vs. time dif- 
— u; in the ID FA model. Top right and bot- 
rescaled response Eq. (|l]) and rescaled autocorre- 
lation function Eq. (|^ as a function of scaled time difference 
{t — w)e~^^'^. Bottom right: out-of-equilibrium FD relations 
in the ID and 2D FA models obtained from parametric plots 
of response vs. difference of connected autocorrelation. 

Having expressed the out-of-equilibrium response and 
correlation functions in terms of the equilibrium ones, 
we can now find the FD relations. Since we know that 
the equilibrium functions satisfy the FD theorem (FDT) 
1, Txeq(r) = Ceq(O) - Ceq(T), from Eqs.(|) and {§ we 
obtain, 



(3) 



where Cc{t, t') = C{t,f) — c{t)c{t') is the connected auto- 
correlation function |25|. This means that for all low T 
and w > 1 there is a unique dynamical FD relation which 
corresponds to the out-of-equilibrium response and au- 
tocorrelation obeying FDT. We illustrate this for vari- 
ous temperatures and waiting times in Fig. 1 (bottom 
right) for the ID and 2D models. Notice that this re- 
sult does not mean that the system is in equilibrium: 
one-time quantities, like the defect concentration c(i), 
are orders of magnitude away from their equilibrium 
values; and, moreover, from the plateau onwards, the 
system develops a strong nearest-neighbour repulsion, 
{dni{t)dni±i{t)) ~ —c(t)[c{t) — Ccq], which is a manifesta- 
tion of the fact that the dynamics explores mostly local 
minima, and eventually vanishes in equilibrium. In fact, 
the response to a uniform field and the autocorrelation 
of the magnetization are not related by any sensible FD 
relation. 

Let us turn now to the more interesting case of the two- 
dimensional square-plaquette model, which consists of a 
system of Ising spins, cr^ = ±1 (i = 1, . . . , N), in a square 
lattice, with ferromagnetic interactions between quartets 
of neighbouring spins in the vertices of the plaquettes of 
the lattice, H — ~J2ijki£a ^i'^j'^k<7i- This model is a 
special case of the eight vertex model with trivial 
thermodynamical properties, but whose single spin-flip 
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FIG. 2. Top left: concentration of defects c{t) as a func- 
tion of time after the quench to T in the 2D square-plaquette 
modeL Top right: integrated response w) as a function of 
time difference t — w. Bottom left: rescaled form of response 
for early times, Eq.(^. Bottom right: rescaled response for 
late times, Eq.(j^, vs. rescaled time c{t){t — u})e~'^^'^ . 

dynamics is glassy [^,^. The model has a dual repre- 
sentation in terms of noninteracting plaquette variables, 
hi = atajakO'i- Since reversing any spin corresponds to 
inverting the four plaquettes to which it belongs, the 
spin-flip dynamics becomes a constrained dynamics in 
the plaquette representation |^ . 

Isolated excitations or defects (hi = —1) in the pla- 
quette model are stable, and have to overcome an energy 
barrier to move. In contrast to the FA model, however, 
pairs of neighbouring defects can move at no energy cost: 
a horizontal (vertical) pair can diffuse freely in a verti- 
cal (horizontal) direction, while a diagonal pair cannot 
diffuse but is free to oscillate. These excitations play an 
important role. An isolated defect moves by creating a 
diffusing pair, so that AE = 2. The single activation bar- 
rier means that this model also behaves as a strong glass. 
Fig. 2 (top left) displays the decay of the concentration 
of defects c{t) = N'^ J2iMt)), where n, = (1 - bi)/2, 
as a function of time after a quench from a random state 
to low temperatures ||2|] . As a consequence of the pres- 
ence of the diffusing pairs the plateau is reached more 
slowly than in the FA model. The activated regime cor- 
responds to times t > 10^ and T < 0.25 in this case. At 
the plateau the defect concentration is cp = 0.089, that 
of moving pairs is 0{e~^^'^ cp), and that of defects which 
belong to oscillating pairs is finite, dp = 0.026. 

In analogy with the FA model, we now consider the lin- 
ear response of the excitations n.i , instead of that of the 
spins ai which in this class of models do not obey any 
systematic out-of-equilibrium FD relations |l|]. In Fig. 
2 (top right) we show the integrated response function of 
individual excitations, x(t, w), for w > 10^ and low tem- 
peratures. The response is again non-monotonic, but in 
this case it also presents an intermediate saturation at 
early time lag. This is a consequence of the existence 
of two well separated relevant timescales. The early be- 



haviour is due to the fast response of the oscillating pairs, 
which exist in finite number in the activated regime. This 
means that we expect, to leading order in e^^/"^, the scal- 
ing form of the response function for low T, w > 10^, and 
t- w ~ 0(1) to be. 



Txit,w)^d{t)fit~w), 



(4) 



where d{t) is the number of excitations which belong 
to an oscillating pair as a function of time, and /(t) is 
an increasing, temperature independent, function, which 
should saturate at a half, /(oo) = 1/2, given that each 
oscillator can be in two states. The collapse under (^) of 
the response curves for early time lags is given in Fig. 2 
(bottom left). After the saturation of the oscillators, the 
response of the isolated defects takes over, so we expect 
the response function to be proportional to c(t), and the 
time lag to scale with the relaxation time, so that for 
t- w > 1, 



Tx{t,w) ^-d{t) + c{t)g 



c{t) 



p2/T 



(5) 



where g{T) is a temperature independent function with 
5(00) = 1, and the relaxation time r^ez ~ e^/"^/c(i) for 
the diffusion of isolated defects scales with time through 
c{t) to account for the necessary number of diffusing pairs 
created to reach another isolated defect. This rather un- 
expected behaviour comes from the restricted ID mo- 
tion of the diffusing pairs and leads to a relaxation time 
Trei ^ e'^/"^ in equilibrium. In Fig. 2 (bottom right) we 
show that (^ accurately scales all the response curves for 
long time lags. Similar scaling relations can be obtained 
for the autocorrelation function. 

The existence of two relevant timescales for the out- 
of-equilibrium dynamics naturally leads to the question 
of whether response and correlation obey nontrivial FD 
relations. In Fig. 3 we present the FD plot obtained para- 
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FIG. 3. Out-of-equilibrium FD relations in the 2D 
square-plaquette model. The inset gives the FD plots for 
various T and w. The main figure is the average of all the 
data. The error bars are comparable to the size of the sym- 
bols. A linear fit to the second part of the curve gives a slope 
X = 0.81 and breaking point Y, = 0.013. 
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metrically from the integrated response w) against 
the difference of the connected autocorrelation Y{t, w) = 
Cc{t, t) — Celt, w) for different temperatures and waiting 
times. The remarkable thing is that the FD relations for 
all the values of T and w fall in a unique master FD plot. 
The FD curve has the piecewise linear structure char- 
acteristic of discontinuous mean-field models Q , with a 
single breaking point = 0.013 and slope x = 0.81. 

Several comments are in order. The autocorrelation 
difference Y{t^ w) is playing the role of (one minus) the 
overlap between configurations in this physical aging situ- 
ation. The dependence of the response on time is through 
this difference, x(*j'"') = x[^(*i Moreover, while 
Y(t, w) is a non-monotonic function of t, it is a mono- 
tonic function of w, so the FD relation can be rewritten 
TxiY) = lo X{y)dy, where X{y) = 1 for y < K, and 
X{y) = X for y > F*. Notice also that for fixed t and 
varying w, the range of Y is determined by c(i), i.e., 
Y{t,w) G [0, c(t) — c^(i)]i ^iid given that for all temper- 
atures under consideration Yf,t^{oo) < Y^, all equilibrium 
curves are contained in the master plot of Fig. 3. 

The piecewise linear FD plot is purely a dynamical 
effect — the static properties of the model are trivial. This 
means that violation of FDT does not necessarily entail 
an RSB-like phase transition as in disordered mean-field 
models 0. The configurations at the plateau seem to 
determine the structure of the FD curve. Each of these 
is composed of cp excitations, dp of which belong to an 
oscillating pair. The breaking point Y^ = dp/ 2 is given 
precisely by the oscillation of all of these pairs. Higher 
values of Y correspond to pairs of configurations mutu- 
ally accessible only through activation. It is harder to 
understand the slope x. The fact that it is independent 
of T means that it is not directly related to the entropy 
of configurations of isolated defects at the plateau, which 
can be calculated from the hard-square model ||2^ . 

The results of this work suggest that nontrivial FD re- 
lations may also exist in the activated regime of more 
realistic strong glasses than the ones studied here, and 
in the physical aging regime of other systems with sim- 
ple statics. While this would only be true for an ap- 
propriately defined class of observables, consistent with 
the remarks of for this class FD relations would be 
well defined and unique. For the simple models consid- 
ered here it corresponds to observables constructed out 
of local energy excitations, which are 'orthogonal' to the 
ones responsible for the activated relaxation, the defect 
concentration in the case of these systems. It would be 
interesting to understand the FD ratio X(y) in terms of a 
geometrical, rather than equilibrium, probability distri- 
bution for the dynamical configurations, in line with the 
ideas of 'Edward's measures' put forward in [|l^. And it 
would also be important to extend these results to the 
more difficult case of fragile systems. 
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